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Abstract 

Topological insulators can be characterized alternatively in terms of bulk or edge prop- 
erties. We prove the equivalence between the two descriptions for two-dimensional solids in 
the single-particle picture. We give a new formulation of the Z2-invariant, which allows for a 
bulk index not relying on a (two-dimensional) Brillouin zone. When available though, that 
index is shown to agree with known formulations. The method also applies to integer quan- 
r~| tum Hall systems. We discuss a further variant of the correspondence, based on scattering 

theory. 

a 

' ' 1 Introduction 

^ Topological insulators are a type of materials that behave as ordinary insulators in the bulk, in 

that they exhibit an excitation gap, whereas the edge has robust, gapless modes. They have 
been theoretically predicted as a class fT5] and as a specific compound [2j , in which the efi'ect was 
then observed [181 HI]. I^i analogy with quantum Hall (QH) systems, the presence of edge states 
is a robust property of the system, in that it is stable under moderate changes of parameters. 
That calls for an explanation in terms of topological invariants. 

As pointed out in [T3], a key feature of topological insulators is fermionic time-reversal 
symmetry. It was shown that two-dimensional time-reversal symmetric insulators admit two 
topologically distinct phases: The phases can not be deformed into one another as long as the 
bulk gap and the symmetry are there. In one phase the insulator is an ordinary one; that is, it 
does not carry currents in the bulk nor at the edges. In the other phase edge states are present. 
They come in pairs and their currents carry opposite charge and spin. This phase is called the 
quantum spin Hall (QSH) phase and its signature is a nonzero spin current, which in contrast 
to the QH current is not quantized, as a rule. For integer quantum Hall systems the topological 
invariant is a Chern number, and can be read off from the value of the Hall conductivity; for 
quantum spin Hall systems the topological invariant is not expressing the value of the spin 
current, but the parity of the number of edge states. 

In this paper we prove the bulk-edge correspondence for two-dimensional topological insula- 
tors and for independent particles. For short we introduce a Z2 bulk topological invariant (to be 
shown equivalent to others), and we show that it is equal to the number of pairs of edge states 
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modulo 2. For sure, that duality has been discussed in the literature, see e.g. [8], but we main- 
tain that there is room for a strict mathematical approach; just as for quantum Hall systems, 
where Laughlin's argument has gained in precision and detail by the subsequent mathematical 
discussion. We should though mention |23|[T]. whose results are further compared to ours below. 

The paper is organized as follows. In Section 2.1 we consider a general class of two- 
dimensional single-particle lattice Hamiltonians, with and without edge. They are symmetric 
under fermionic time reversal and hence potentially describe topological insulators. One note- 
worthy feature is that periodicity is postulated only in the direction parallel to the edge. In 
Section 2.2 we formulate the duality in its most basic version, in a precise though still prelim- 
inary form. The purpose of Section [3] is to show how graphene and related model materials, 
which have been considered as candidates for topological insulators, fit into our scheme. As a 
side remark, we give a simple new proof of the absence of edge states for armchair boundary 
conditions. In Section [4] we discuss some abstract vector bundles on the 2-torus. By considering 
their sections and transition matrices, we classify them in terms of a Z2-invariant, the formula- 
tion of which might be new. So equipped, we define the bulk index in Section [5] and state in full 
the basic version of duality (Theorem 5.4). It is immediately followed by the main steps of the 
proof, while technical details are postponed to main part of Section [7j The part of the article 
so far described provides a full and self-contained account of basic bulk-edge duality. It is the 
most general part, and yet makes up for less than half of its length. 

In Section 5.1 we formulate a bulk index for the specific and more familiar case where the 
lattice Hamiltonian is doubly periodic. As it is to be expected, the relevant 2-torus is now 
the Brillouin zone. We discuss how that index arises from the general one, leaving details to 
Section [Sj All the results obtained up to that point have a counterpart in the case of QH 
systems, which we present in Section [6] for illustration, because they are simpler in that context. 
As a matter of fact, for such systems we include an independent, alternate version of the duality 
based on scattering theory and more precisely on Levinson's theorem (Section 6.2), as well as a 
comparison between the two versions (Section 6.3). Related details are found in Section [oj We 
conjecture an analogous alternate version for QSH systems. Finally, Section 7.1 contains some 
results about indices, including a comparison of ours with some of those found in the literature. 

Let us comment on the relation to other work. In |15j an edge and a bulk characterization of a 
topological insulator is given, but without proof of a mathematical link. In [23] a correspondence 
between bulk and edge description is given. There however the bulk is put in correspondence 
with twisted boundary conditions allowing tunneling between a pair of edges; they include open 
boundaries as a special case. Moreover, the bulk invariant used there is physically different from 
a Z2-classification and may not reflect itself in open boundaries. In [1] the bulk-edge duality for 
topological insulators is derived in the same sense as ours, among other results. The setting, 
however, is not quite as general as ours, as far as the side of the bulk is concerned. As the 
edge is concerned, an index related to the Maslov index is defined and shown to be equivalent 
to that of fl5]. As for the bulk, it is assumed that the Hamiltonian is the perturbation of one 
commuting with spin, and hence consisting to several copies of a quantum Hall insulators. It 
is assumed that the gap remains open as the perturbation is switched on. Thus spin resolved 
Chern numbers can be defined by homotopy [22\ 177] and the duality inferred from that of the 
quantum Hall case [26]. 

Further, assorted comments are: Some of our examples are discussed in [Ij in a similar vein. 
Related indices for topological insulators are discussed in [H [8l [HI [l5l dU [25] . A similar choice 
of a torus and bundle as ours in Sect. 2.2 is found in [28]. The results of Section 5.1 depend 
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on the analytic properties of band functions [T7] . In [T2] bulk-edge duality for Hall insulators 
is pinpointed at the birth of edge states of of band edges. This insight is the reason for using 
Levinson's theorem, though the method is otherwise different and the result more general. 
Last but not least, the existence of a complementary approach to gapped systems, including 
interacting ones, should be mentioned. It is based on effective, topological field theories as a 
tool to explore the response of a system in the limit of low frequencies and long wavelengths; 
see [221 H El EOl El E] for early examples. 

2 Setting and results 

We shall introduce a class of bulk, resp. edge, single-particle discrete Schrodinger operators, 
by which we describe insulators, topological or otherwise, which extend over a plane, resp. a 
half-plane. We conclude the section with a loose description of the results. 

2.1 The Schrodinger operators 

Consider, at first and in less than final generality, a tight binding Hamiltonian with nearest 
neighbor hopping on the lattice Z x Z, resp. N x Z. The Hamiltonian is assumed to be periodic 
in the direction along the edge, but not necessarily across it. The period may be taken to be 
equal to one without loss: In fact sites within a period may be regarded as labels of internal 
degrees of freedom, among others like e.g. spin. We may thus perform a Bloch decomposition 
with respect to the longitudinal quasi- momentum k e := M/27rZ, which remains a good 
quantum number even in presence of the boundary. 

We so end up with a family of Hamiltonians H[k) defined on the one-dimensional lattice 
Z 3 n, resp. N, and acting on wave-functions ■i/'n £ C^, where N is the number of internal 
degrees of freedom. 

More generally, these objects are stated as follows. 

Definition 2.1. [Bulk Hamiltonian] The Hamiltonian, acting on ip e C"^ {Tj; C'^ ) and parametrized 
by k e , is 

{H{k)ll:)^ = A{k)iJn-l + A{kY^n+l+Vn{k)^n, (n G Z, £ C^) • (1) 

The potential Vn{k) and the hopping matrices A{k) are N x N matrices having a -dependence 
on k, uniformly in n. We assume Vn{k) = Vn{k)* , where * denotes the matrix-adjoint, and that 
A{k) e GL(A^). (Recall that GL(A'^) M7v(C) consists of invertible matrices of order N.) 

We then consider the restriction of the Hamiltonian toN = {l,2,...} (we find it convenient 
to omit zero), while allowing for changes within a finite distance ng ^ from the edge. 

Definition 2.2. [Edge Hamiltonian] The Hamiltonian, acting on t/j e ^^(N;C^), is 

= A{k)i;n-i + ^(A:)*Vn+i + V^{k)^Pn , (n e N, ^„ e C^) , (2) 

where satisfies the above properties of V, as well as 

V^{k) = Vn{k) , (n > no) . (3) 

Moreover we assume the Dirichlet boundary condition, meaning that for n = 1 Eq. is to be 
read with ipQ = 0. 
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Remark 2.3. The Dirichlet condition is by no means special. Since 

any boundary condition ifjQ = htpi, (A{k) e M]\j{C)) defines a self-adjoint Hamiltonian if 
{AK)* = AK. That amounts to the Dirichlet condition after adding 5niAK to Vn. 

Contrary to what the above motivation might suggest, the Hamiltonians ([l|[2]) are prompted 
by more than just the square lattice. In the next section we will show that several models based 
on the honeycomb lattice, which have been considered [15] in relation with the quantum spin 
Hall effect, fit the scheme. Moreover, one-dimensional spin pumps [8] also match the description, 
with k playing the role of time. 

The topological classification applies to insulators that are invariant under odd (or fermionic) 
time-reversal symmetry. In the sequel we specify these notions. 

Definition 2.4. [Time-reversal symmetry] The symmetry is a map : — > with the 
following properties. 

i) @ is antilinear and = —1; 

ii) e*e = 1; 

Hi) For all k e , 

H{-k) = eH{k)Q-^ , (4) 
where G also denotes the map induced on £^(Z;C^). Likewise for HK 

As a result, N is even. In the models described in the next section, properties (i, ii) arise 
from the time-reversal of a spin-^ particle, and (iii) from the symmetry of the Hamiltonian. 

The Bulk Hamiltonian of an insulator is supposed to have a spectral gap at Fermi energy /i 
and for all k: 

H^a{H{k)), (keS^). (5) 

Typically -ff" does not satisfy the gap condition. In fact, while for the essential spectra we have 
aess{H^{f^)) o'ess{H{k)), the edge Hamiltonian may have discrete eigenvalues crossing /x for 
some values of k. 

2.2 The main result in brief 

We can informally introduce two indices, X, e {+!}, defined in terms of H and ij", respectively. 
On the circle the involution k ^ —k has two fixed points, k = 0, vr. The edge index l'^ = (—1)'^ is 
the parity of the number n of those k e [0, vr] at which an eigenvalue of H'^{k) equals //, at least 
if the eigenvalue crossings are simple. The index X requires more explanation. By definition 
H[k) does not have eigenvalues z ^ a{H{k)); yet we may regard the Schrodinger equation 

{H{k)-z)^ = 

as a second-order difference equation inn e Z. As such, it has 2N linearly independent solutions, 
but we focus attention on those which decay as n — > -l-oo. They form a linear space, £'z,fc, of 
dimension N. By (i, iii) we have 

{H{-k) - z)e^|J = e{H{k) - z)tl; , 
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and thus E^-k = QE^^k- We choose a reflection symmetric complex contour, 7 = 7, encircling 
the part of the spectrum of H{k) lying below /i, and we set T = 7 x S"^. We so have: (a) an 
involution {z,k) >—>■ {z,—k) on the torus T; (b) a vector bundle with base T and fibers E^^k', 
which (c) are compatible with in the stated sense. To any vector bundle with these features, 
including (i), but irrespective of the concrete definition of its fibers, an index will be associated 
in Sect.|4j The bulk index I is that index for the particular bundle arising from H as described. 
The main result of this work is that I = iK 




Figure 1: The torus T is the product of the loop 7 in the z-plane with the fc-axis (both red); it 
is to be glued along the dashed loops. Thick lines (in black) delimit parts of the bulk spectrum 
a{H(k)), as a function of k. Curves (in green) connecting them are discrete eigenvalues £{k) e 
a{H'^{k)), representing edge states. They intersect the torus only at crossing points (dots) along 
the Fermi line {z = fi} x (in blue), one half of it (thick) being relevant for the index XK The 
spectra are symmetric in A; 1-^ —k. 



It should be noted that the result relies on just one translational symmetry. 



3 Examples 

We show how to obtain Hamiltonians [2]) from the Schrodinger operator on the honeycomb 
lattice, which models graphene in the single-particle approximation. We will consider two types 
of boundary conditions: zigzag and armchair. It is known that the spectrum of the Schrodinger 
operator depends on which boundary condition is chosen, [201 [H]- Iii particular, a zigzag 
boundary implies the presence of zero-energy edge states, while they are absent for the armchair 



boundary. In Example 3.1 we give a new proof of the last statement. 



3.1 Graphene 

Let A = u As be the (infinite) honeycomb lattice and its bipartite decomposition into two 
triangular lattices A^, A^. Upon fixing an origin, they are 

A^ := {n = niai + ^202 | (ni,n2) £ Z^} , Ab := Aa + d , 

IT - 1 

|ai| = |a2| , Z(ai,a2) = - , S = -{ai + 02) ■ 
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Any site in A a has three nearest neighbors in A^, shifted by 5 or by another equivalent vector. 
See Fig. [2j The model for graphene is simply the Schrodinger operator Hq for a particle hopping 
between nearest neighbors (with hopping parameter —t). 




Figure 2: The honeycomb lattice. Coordinates of lattice sites are relative to the basis {ai,a2}. 

The (Bravais) lattice of translations of A is A^. We reduce A to A^^ by dimerizing neighbors 
shifted by 5. We retain the position n of the ^-site as that of the dimer, and the values of the 
wave function at the two sites as pseudospin components: 

■■= (^l) e . (6) 
With these notations, the Schrodinger operator Hq takes the form 

\V^ni,n2 + l ~ V^ni + l,n2 ^ni,n2/ 

The expression for Hq could of course have ended up differently. Let us name the choices 
underlying its construction: a sublattice shift vector (above: 6), defining the dimer, and two 
primitive lattice vectors (ai and 02), defining adjacency between dimers. For instance for = 
ai, a'2 = ai + 0,2, and hence n\ = ni — 77-2, n'2 = ^2, Eq. ^ becomes 

('^f„,/,^ _ i ( '^ni,n2 + V'ni + l,n2-l + '^ni-l,n2\ 

l^O'i/'j - -t I , 4 I A . l A )> 

after dropping primes. See Fig. |3j 

The Bloch decomposition w.r.t. n2 e Z 

r fiu 

^n= 7^e^'"^^„i(A:), (V'„i(A;)eC2) (8) 
fibers the Hamiltonian, (//o'0)ni(^) = -f^o(^)V'ni (^)- feet we obtain from Eq. ([T]) 

{Ho{k)^)^ = AQ{k)lPn-l + AQ{kyiPn+l + VQ{k)%bn , 

where we set ip = il^{k), n = ni and 
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Figure 3: The honeycomb lattice with zigzag a) and armchair b) boundary conditions. The 
lattice has been rotated, so as to have ni fixed along the vertical. 

Likewise for HQ{k): 

{Ho{k)^)^ = Mk)i'n-i + MkTi'n+i + Vo{k)i;n , (10) 

where now 

Mk) = ~t J) , Voik) = (J J) . 

It should be noted that only the second Hamiltonian satisfies the condition A{k) e GL(A^) in 
Eq. ([T]). The two equivalent bulk Hamiltonians no longer are once they are turned into edge 
Hamiltonians by means of the Dirichlet boundary condition. Indeed, they correspond to (a) 
zigzag and (b) armchair boundary conditions, respectively. See again Fig. |3j 

Proposition 3.1. (i) The Hamiltonian Hq has E = as an eigenvalue. Actually, -ffo(^) 

for k e (27r/3, 47r/3). (ii) The Hamiltonian Hq does not have any eigenvalue, i.e. its spectrum 

is purely continuous. 

The result is known [9l [20] , but perhaps not the argument below for (ii) . 

Proof, i) Related to AQ{k) being singular, HQ{k)ip = reduces to the first order equations 

+ (1 + e-*'=)V'^ = , + (1 + e''')^^J^ = , (n e N) . 

The boundary condition ipo = implies = 0, but is ineffective for ip^, for which there is a 
non-trivial solution in £^(N), as long as |1 + e*'^! < 1, i.e. for k in the stated range. 

ii) Every second dimer of the armchair boundary is exposed (see Fig. [3] right) and the axis 
containing them is ni = 0. Let us consider, on the whole lattice A, the reflection n i—> rn of 
lattice sites about the axis ni = in Fig.jsj as well as parity, {P'iIj).^ = iprn- Clearly, [P, i/o] = 0. 
Any odd eigenfunction of Hq satisfies the Dirichlet boundary condition on the line ni = 0, 
whence its restriction ^ to ni > Q defines one for H^. What matters more is that the converse 
is true as well: The odd extension tp ol ip satisfies the Schrodinger equation on ni = 0, and 
hence at all ft. The conclusion is by noting that Hq has no eigenvalues; actually it has purely 
absolutely continuous spectrum. This is known and by the way has a short proof: It suffices to 
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establish the property for the fiber Hamiltonian Ho^k). Since the latter has real analytic Bloch 
eigenvalues A(g), absolute continuity holds unless A(g) is constant in q. That, finally, is ruled 



out by Ao{k) e GL(A^) in Eq. (10) and the remark below. 
Remark 3.2. (l^, Prop. 7) Suppose the Hamiltonian on £^(Z;C^) 



n 



(note: V independent of n) has a Bloch eigenvalue \{q) independent of the quasi-momentum q. 
Then A is singular. 

Proof Let A be that eigenvalue. Then det(yl2;^^ + A*z + V - X) = foi z = e**?, {q e M) and, 
by analyticity, for z ^ 0. Letting z — > gives det A = 0. ■. 

3.2 The Kane-Mele model 

In the next example we discuss a model with odd time-reversal symmetry. Wave functions are 
now of the form ip = {ipri)nei?i where n labels dimers, see Fig. @ and V'neC^^C^. The spin of 
the particle is represented by the Pauli matrices ai , (T2 , as acting on the second factor . The 
time-reversal operator is = — i(l ® o"2)C, where C denotes complex conjugation, which is the 
usual expression for a spin-j particle. 

Let Hq be the Schrodinger operator on the honeycomb lattice defined in Eq. ([7| . We consider 
the operator 

Hkm = Ho®1 + Hi^as , 

where the last term implements a spin-orbit coupling. There Hi describes hopping between 
next-to- nearest neighbors; it thus acts diagonally on the pseudospin components of Eq. and 
specifically as 

with 

rti — l,n2 + l "I" ^ni + l,n2 ~^ V'ni,n2 — l) ^(V'ni,7i2 + 1 ~^ ^ni + l,n2 — 1 ~^ V'rai — 1,712)) ' 

The grouping of terms reflects that, starting from a given ^-site, its next-to- nearest neighbors 
are reached by turning right or left at a i?-site. The turns go along with phases +i = e-*'^/^ 
modeling a magnetic flux tt/2 through any (positively oriented) triangle ABA; yet the total flux 
through the hexagonal cell vanishes, since no phases are associated with the bonds forming its 
boundary, see Eq. ([T]). Likewise for B and A interchanged, resulting in the sign in (11). Indeed, 



for a given dimer, the next-to-nearest neighbors form again dimers, but the turns linking A, 
resp. S-sites are opposite. 

The model is the special case of the Kane-Mele model, |15j . where two further terms (a 
Rashba term and a staggered chemical potential) have been set to zero. Even so, it exhibits a 
non-trivial topological phase, i.e. X = —1, for small t! jt > and fi = 0. In fact, it is the direct 
sum of (a special case of) the Haldane model [TO], Hn = Hq + Hi, and of its time-reversed copy 
Hq — Hi. In such a situation we have I = (— l)'^, where J\f is the integer associated with the 
quantum Hall effect of Hn. In the stated regime, a band gap of order 0(|f'|) opens, and M = 1, 

m. 
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We conclude the example by giving the fibers of //km w.r.t. the Bloch decomposition 
It will suffice to do so for Hn: Instead of Eq. Q we have 

4 Time-reversal invariant bundles 



The purpose of this section is to define the index of bundles of the kind mentioned in Sect. 2.2 



see Def. 4.8 below, as well as some auxiliary indices; and to formulate some of their properties. 



We refer the reader to Sect. [7] for the proofs of the lemmas stated here. 
4.1 The index of endpoint degenerate families 

We shall define an index for certain families of points on the unit circle. We will give the most 
general definition, in order to make evident the stability of the index under homotopy. At the 
end of this section we give a procedure to compute it in a more restrictive setting, which is 
suitable for our application. 

Following ([IS], Sect. II. 5. 2) we consider unordered A^-tuples Z = {zi , . . . , zjy) of complex 
numbers Zi e C A distance is defined by 

d{Z' , Z) = minmax \zn — z'^l , 

n 

where the minimum is taken over all possible relabellings of Z or Z' . Let us recall (|16j. 
Thm. II. 5. 2): Given a family Z{x), (x e [a, ^]), which is continuous w.r.t. d, there exists a 
(non-unique) continuous labeling Z{x) = {zi{x) , . . . , zn{x)). 

In the following we will consider A^-tuples with Zi e S^. Any continuous labeling Zi{x) 
induces continuous arguments Oi{x), i.e. Zi{x) = e^^i^^) . Two simple observations are in order: 
(i) ^Z):=(2^)-il]f^ -i is independent of the choice of arguments, as well as of labeling. 

In fact ]^^]^^;j(a;) is independent of the latter, (ii) Any family Z{x) with Z{a) = Z{b) has a 
winding number given as 

J\f{Z) := w{Z) . (12) 

We next consider endpoint degenerate families Z{-): for x = a,b each z e Z[x) occurs with 
even multiplicity. We may concatenate such a family with one, Z(x), (x e [b, c]), such that 
Z{b) = Z{b) and Z[c) = Z{a), while keeping Zi{x) even degenerate. Clearly, Z is not unique. 



but by (12) 



w{Zi) -w{Z2) e 2-Z . 

By the same reason, M{Z#Z) is an integer; by M^Z^j^Z) = w{Z) + w{Z) it is determined 
mod 2 by Z. 

Definition 4.1. [Index] We set 

X{Z) = ^ (13) 

as the index of endpoint degenerate families. 
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Consider now the special situation where there is z e such that z e Z{x) occurs only at 
finitely many which moreover are simple crossings: z = Zj{x) for a single j, and z'-{x) 7^ 0. 
Then 

X{Z) = {-IT , (14) 

where n is the number of crossings of z. Indeed, one can choose Z without crossings of z. Then 
N{Z^Z) is the number of signed crossings of z by Z{-) but that qualification is irrelevant for 
parity. 

4.2 The index of Kramers families of matrices 

Let £ be the matrix of even order N given by the block diagonal matrix with blocks 

'0 -1' 



■} 

C : the complex conjugation, and Go = eC the standard time-reversal on C^. 

Suppose T e GL(A^) satisfies 

BoT = T^^Qo . (15) 

Then the eigenvalues of T come in pairs A, A^^ with equal algebraic multiplicity; in particular 
their phases z = A/|A| are even degenerate. Indeed, @o{T — A)" = T^"(l — AT)"0o, as seen 
inductively for n = 0, 1, . . .. 



Definition 4.2. [Kramers property] We call Eq. (15) the Kramers property. We say a family 
T{ip) e GL(A^), which is continuous in ^ f ^ n, has that property if the endpoints T(0) and 
T{tt) have it. 

The repeated eigenvalues \i{'-p) of T{(p) form a continuous family in the sense of the previous 
section ([16|, Thm. II. 5.1), and so do the Zi = Ai/|Ai|. Moreover, Z[ip) = [zi^^p),. . .Z]y{ip)) is an 
endpoint degenerate family. 

Definition 4.3. [Index] We set 

I{T) = I{Z) (16) 

as the index of a Kramers family T. See Eq. (13). 

All it in fact takes for the definition is the endpoint degeneracy of Z and not the stronger 
Kramers property of T. However we shall not need such an extension. 

Lemma 4.4. Suppose 

T2{^) = M^{ip)Ti{^)M+{ip)-' 

with continuous M+{(p) e GL(A^) (0 ^ ip ^ tt), as well as 

M_(0)eo = eoM+(0) , M_(^)eo = eoM+(7r) . (17) 

Then M{(p) := M_(c/?)M+((/3)^-'^ has the Kramers property, and T2{ip) has it iff Ti{{p) does. If 
so, 

X{T2) = Z{Ti)Z{M) . (18) 



The claims are of immediate verification, except for Eq. (18). However in the special case 
that M+ ((/?) = M_ {if) the equality 

X{T2)=X{Ti) (19) 

is also immediate, because Ti{ip), T2((/?) then have the same eigenvalues. That case suffices for 
the basic result on bulk-edge duality. 
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4.3 The index of time-reversal invariant bundles 

Let = M/27rZ be the circle and T = x S'^ 3 tp = {ipi,ip2) the torus with involution 
T : 99 H-> —{p. It has four fixed points: ipo = (0,0), (0, tt), (7r,0), (7r,7r). Let E he a complex 
vector bundle with base T 3 99 and fibers of dimension N. We say that E is time-reversal 
invariant if there is a map Q : E —>■ E with = —1 and 6 : E^ >—>■ E^^ antilinear. 

We also consider the associated frame bundle F{E) with the following operations induced 
on frames v = {v\,. . . ,vn) e P{E)(p: 

• right multiplication by M e GL(Ar): 

M : F{E)^ ^ F{E)^, v ^ vM , (20) 

with {vM)j = Xiili ViMij. Any two frames are so related by a unique M. 

• e : F{E)^ F{E)r^, V ^ Qv, with {Qv)i = Qvi. Note that 

Q{vM) = {&v)M . (21) 

In order to classify time-reversal invariant bundles we consider the torus cut along the circle 
{ipi I (^1 = TT = — tt} X S^, ov more precisely T = [— 7r,7r] x S^. Note that only (p = (0,0), (0,7r) 
remain among the fixed points of r, and that any bundle on T naturally defines one on T. 

Lemma 4.5. [Existence of time-reversal invariant sections] On the cut torus T, there are 
(smooth) sections v : T — » F{E) of the frame bundle (whence v{(p) e F{E)^) satisfying 

v{t(p) = ev{^p)e . (22) 

As a result, they are compatible with right multiplication by M{ip) e GL(A'') iff 

QqM{t^) = M{ip)eo . (23) 

Let v+{ip2) ■= v{±Tr,ip2) be the boundary values of v{(p) along the two sides of the cut, and 
T{ip2) e GL{N) the transition matrix, 

V+{ip2) = v-{^2)n^2) , {^2 £ S^) . (24) 
Lemma 4.6. [Time-reversal symmetry of the transition functions] 

Go 'r(-(^2)eor(<^2) = i . 

In particular, T{ip2) has the Kramers property on the interval ^ 992 ^ tt, which hints at a 
the possibility of assigning an index to the bundle. 

Lemma 4.7. [Independence of the index from the section] Let v^'^\ (i = 1,2) be time-reversal 
invariant sections on the cut torus T, and let Ti{(f2) the corresponding transitions matrices 
across the cut. Then 

I{Ti) = X{T2) . 
We may thus proceed to the following definition. 
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Definition 4.8. [Index] We set 

I{E) = I{T) (25) 



as the index of a time-reversal invariant vector bundle E over T. See Eq. (16). 



Such bundles can hence be distinguished in two topologically distinct classes, according to 
the value of the index. 

We also retain the following remark, which appears as a byproduct of Lemma |4.5[ 



Remark 4.9. Let Q : satisfy the time-reversal symmetry condition (i) of Def. 2.J^ 

Then there is a basis v = [vi, . . . , vn) ofC^ such that v = Qve. If condition (ii) applies too, the 
basis can be chosen orthonormal. Letting U map v to the standard basis of , a restatement 
is UQ = —QqU with U unitary. 

5 The bulk-edge correspondence 

It pays to look at first at wave-functions 'i/' = (V'n)neZ just sequences, i.e. ^ e C := Z x C". It 
can be shown that for any z e p{H{k)) in the resolvent set, the Schrodinger equation H{k)^ = zip 
has N linearly independent solutions which are square- integrable at n — > -l-oo. Let E^ k ^ C the 



linear space they form. As explained in Sect. 2.2 we obtain a vector bundle 



E = {{{z,k),il;)eTxC\ijeE,^k} (26) 

over the torus T = ^ x S^, which enjoys the property of being time-reversal invariant. The 
following definition is thus natural. 

Definition 5.1. We set 

I = I{E) (27) 



as the bulk index. See Eq. (25). 



The edge index has been loosely introduced in Sect. 2.2 In preparation for precise definition 
let us take a closed interval I containing the Fermi energy /i in its interior and such that 
a{H{k)) n I = for k s S^. The spectrum of H'^{k) is discrete and finite in I (uniformly in k), 
and the eigenvalues are locally given by branches £i{k) which are in k ([H], Thm. II. 6. 8). 
By possibly adjusting fi we can arrange that if £i{k:^) = fj, for some /c* then £'j(/c*) 0. We note 
that the number of crossings at A;,^ = or A;,^ = vr is even by Q . 

Definition 5.2. Let n be the number of crossings k^ e [0,7r], with those at endpoints counted 
half. Set 

Xt = (-1)" 

as the edge index. 

The next remark about the edge index is inessential for the following main result. 

Remark 5.3. The index X" could be also defined in terms of a Kramers family. To this end 
let /(e) be a continuous real function with f{e) =0, (e < I) and /(e) = 1, (e > I), where 
I is the aforementioned interval. Then, in view of Qg{H)@^^ = g{QH&^^) , the operator 
T{k) = exp(27ri/(i/''(A;))) has the Kramers property, albeit w.r.t. Q, which is however irrelevant. 
The definition l'^ = T{T) would be legitimate, since that ofw{Z) in Sect. 4-1 extends to countable 



families Z of points Zi{x) G , as long as only finitely many (uniformly in x) are ^ 1. It would 



agree with Def. 5.2 by Eq. (I4) 
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Theorem 5.4. [Bulk-edge correspondence for topological insulators] Let the Bulk and the Edge 
Hamiltonian be as in Defs. 2.1 and 2.2. Assume the time-reversal symmetry conditions of 
Def. 2.4 and the gap condition Then 

i = iK 

We can give the main steps of the proof right away. For given [z, k) e p{H{k)) x S"^ we 
consider, besides of E^^ki also the linear space E"^ ^ of solutions ■0" = (V'n)neZ decaying at n ^ cx) 
of the Schrodinger equation H'^{k)il)^ = ztp'^, without imposing boundary conditions at n = 0. 
Eq. ([3]) establishes a bijection 

E^,k - Ef,fc , ^ - V'" (28) 

determined by ^pt, = ifjn for n > uq. 

Frames ^ e F{E)z^k consists of A^-tuples ^ = (-01, ■ ■ • , V'iv) of linearly independent solutions 
ipi e Ez^k- Note that the index i does not denote the lattice site n e Tj. Since £ C^, we 
may equivalently say: ^' = (^'n)nez with e M]\f{C) belongs to F{E)z,k iff ^' is a solution of 
H{k)^ = z'^ decaying at n ^ +00 , which is moreover fundamental in the sense that for any n 

^'„a = 0, ^„+ia = ^ a = , (a e C^). (29) 



The bijection (28) induces one between frame bundles, F{E)z^k — >■ F{E'^)z^k-, ^ ^ It is 
manifestly compatible with the right action (23) of GL(A^). The next lemma rests on the 
bijection. 



* 



Lemma 5.5. i) A point {zn:,kit:) e T has det ^'q = for some (and hence all) e F{E)z^^k 
iff z^ e a{H'^{ki^)). If so, then z^ = ^. For a dense set of Hamiltonians near the given 
one, the points k^^ are isolated in and for each of them there is a simple eigenvalue 
branch e{k) with e{k^) = fi, e'{k^) 7^ 0; moreover, 

det ^5^0. (30) 

Density is meant with respect to the topology of the class of Hamiltonians specified at the 
beginning of Sect. \2.1\ 



ii) Let "^{z, k) e F{E)z^k be a section defined in a neighborhood in C x ^ T of any of the 
crossing points {z^ = fx, k^). The family of matrices 

L{z,k) = -¥^*{z,k)A{k)¥Q{z,k) 

has the reflection property L(z, k) = L(z, k)* . Its eigenvalues are thus real for real z. There 
is a single eigenvalue branch l{z, k) of L{z, k) vanishing to first order at (//, k^f): There the 
derivatives dl/dz and dl/dk are real and non-zero. 



Hi) At any of the points (/i, k^:) we have 

dl 



< . (31) 
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iv) As k increases past k^ the eigenvalue £{k) crosses fi as an increasing function if 



dk 



(z=/i,A;=A;*) 



>0, 



and as a decreasing one in the opposite case. 



The lemma, which is proven in Sect. [7j allows to complete the proof of the main result. 
As a matter of fact, only items (i, ii) matter to that end. For later use we mention that no 
use of time-reversal symmetry has been made in the lemma, in the sense that the statement 
about density of the Hamiltonians in (i) holds in either class, with and without that specification. 



Proof of Thru. 5.4- As a preliminary we recall Remark 4.9 At the price of conjugating the 
Hamiltonians by U in the internal space C^, which preserves the assumptions, we may assume 



e = -e, 



• 



(32) 



By density it will suffice to prove the theorem for Hamiltonians as specified in part (i) of the 
lemma. By the first sentence there, we can define a section (z, k) ^ F{E)z^k away from crossing 
points by requiring 

*»(z,A:) = l. (33) 

Clearly, the requirement can not be imposed at such points. In a small reflection symmetric 
disk D = D <^ T containing the generic crossing point {fj,,k^), Eq. (30) allows us to make an 
alternate choice ^{z,k) by requiring 



¥,iz,k) 



1 . 



(34) 



Figure 4: The neighborhood in T (see Fig. [T]) of a crossing point (/i, /c*). Two different local 
sections are defined inside (^') and outside (^') of a small disk D containing (fi,k^). They are 
glued across the two half-circles (^dD)+ by means of matrices M+{k). 

Based on the latter section we shall extend the former to the cut torus T, where the cut is 
the Fermi line {z = fj,} x S^. On the boundary BD both sections are defined, and we have 

^{z,k) = '^{z,k)M+{k) , (35) 



for {z,k) e {dD)+ = dD n {+Imz > 0} and some matrices M+{k) e GL(7V), see Eq. (20), 
which are parametrized by k in the interval J resulting from the intersection of the Fermi line 
with D. The same matrices relate and We may then obtain the announced extension 
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by (re) defining ^{z,k) through Eq. (35) for (z,k) e D+ = D n {+Imz > 01. We observe that 
the section on T so constructed satisfies Eq. (22); in fact its ingredients (33, 34) do by (32) and 
(-eol)e = 1. 

The transition matrix T[k) across the cut differs from the identity only within the intervals 
J. Along such an interval the boundary values of the section are ^'+(A;) = "^{fj,, k)M+{k)^^ , 
resulting in the transition matrix 



T{k) = M-{k)M+{k)-^ , {kel). 



(36) 



We can compute the bulk index T{E) = T{T) using Eq. (14). The theorem then reduces to 
the claim that, for each crossing point, the eigenvalues of T{k) change with k e J from 1 to 1 
without winding, except for a single one with winding number +1. 

It is with the proof of this claim that our choice of sections comes to fruition. In fact, by 



Eq. (33) and the sentence after Eq. ( |35| ) we have 

¥o{z,k) = M+{k) , {{z,k)e{dD) + ) . 

We apply part (ii) of the lemma to the section ^l*, which is as required there. Then L{z,k) 
-A{k)¥Q{z,k) by Eq. (34| and hence 



T{k) = A{k)-^L-{k)L+{k)-^A{k) 



ike I) 



with L+{k) = L{z, k) for [z, k) e {dD) + . Clearly, A{k) can be dropped, as it does not affect the 
eigenvalues under investigation. We have 

L{z, k) = l{z, k)U{z, k) © Z(z, k) , 

where n(z, k) is the rank 1 projection on the eigenspace corresponding to the eigenvalue l{z, k), 
and L{z, k) is regular as a map on the range of the complementary projection. Clearly, l{z, k) = 
0(A) as A := (z — /i, /c — k^) — > 0, but by (ii) above we also have l{z, k)^^ = 0{A~^) on T. We 
use the notation /o and /i for the value of a function / and its gradient {dzf, d^f) at the critical 
point. Then 

L{z, k) = l{z, k)Uo + Lo + Li-A + 0(A2) , 
L{z, k)-^ = l{z, k)-^U{z, k) © L{z, k)-^ = l{z, k)-^{IiQ + Hi • A) + Zq ^ + 0(A) . 

In view of L^{k) = L{z, k) for (z, k) e {dD)+ we compute 



L(z-,A:)L(z,fc)-i = g^no 



ii + (1 - Ho) + 0(A) , 



(37) 



R = l{z, k)-\LoUi • A + Li • AHo) = Z(z, k)-'LoUi • (A - A) = 0(1) . 



In the second line we used -LqIIi +LiIlo = 0, as seen from expanding LII = 0. Using (1— no)ni = 
IIiIIo from n = n^, we find that i? = (1 — no)i?no is strictly block triangular. For small D and 
{z, k) 6 {dD)+ the eigenvalues of ( [37| ) wind as if the error 0(A) is omitted, provided those of the 
explicit part do not vanish. Those in turn remain the same when R is omitted, and are in fact 
equal to 1 with multiplicity A'" — 1, and to l-{k)/l+{k) with multiplicity 1, where l+{k) = l{z, k), 
{{z,k) e {dD)+). The winding number of l-{k)/l+{k) along J equals that of l{z,k) along dD. 
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Since the vanishing at the crossing point is of first order, /(z, k) = li ■ /S. + O(A^), that number 
is 

sgn((5,/)(40) = +1 , 

as claimed. Though the sign is irrelevant (so far), we observe that by (iv, v) it is +1 if the 
eigenvalue crossing is decreasing in fe. ■ 



5.1 The bulk index as an index of Bloch bundles 

A feature of definition of the Bulk index is that it applies to Hamiltonians which are periodic 
just along the edge, see Eqs. ([T 27). If the Hamiltonian is periodic in both directions and hence 
the Brillouin zone two-dimensional, the Bulk index allows for an alternate formulation in terms 
of the bundle of Bloch solutions, as we are about to explain. 

We temporarily suppress the longitudinal quasi-momentum k in Eq. ([T]) and consider Hamil- 
tonians of ^^(Z, C^) of the form 



(38) 



where A e GL(A^) and Vn is periodic in n, i.e. Vn+M = Vn- By once again considering sites n 
within a period as labels of internal degrees of freedom, we may assume that the period is 1. 
This amounts to the replacement of ipn, A, Vn by 




A 



/O 





^\ 




V 



fVo A* 

A ••. 

■•. 

Vo ... 







A 



\ 


A' 
Vm-1 J 



In particular A is singular, unlike A. A Bloch solution {ipn)nei, of quasi-periodicity ^ 0, 

tpn+pM = C^V'n , 

and of energy z is then represented as a solution ^ of 



(39) 



(40) 



(41) 



For K e = ]R/27rZ and ^ = e***, the matrix 'H{C) is hermitian, since now ^ = It thus has 
eigenvalues z = Xi{k) (real and increasingly ordered) and eigenvectors ^'/(k), (/ = 1, . . . ,MN). 
The ranges of the energy curves Xj{K,) are known as energy bands. Let 



Ai{k) = inf A;+i(/t) - sup A;(k) 



(42) 



be the gap between successive bands of H{k), where we temporarily reinstated the dependence 
on A; € S^, implicit in (38). The gap is open if Ai{k) > 0. Let then A; = inffcA/(A;). In 
topological insulators the bands are degenerate for A; = 0, vr, in fact at k = 0, vr. A band thus 
can not be separated from the rest of the spectrum, but a pair of them can. We will assume so 
for the pairs (2j - 1, 2j), (j = 1, • • . , iVo/2), i.e. 



Ao, Aa 



Ano > 



(43) 
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where Nq (even) is the uppermost band below the Fermi energy /i. For that band we actually 
retain the stronger assumption ([s]) of a spectral gap. (It amounts to a positive gap in Eq. (42) 
for I = No even when extremizing jointly in k, n.) 

Definition 5.6. The Bloch bundle Ej of the j-th pair of bands has the Brillouin zone B = 
X 3 {k, k) as base and the span [^2j-i(^) ^)) '^2j{k, k)] cz C^^^ as fibers. 

It should be noted that while the eigenvectors are not smooth in A;, k, their span is (|16]. 
Sect. II. 1.4), since degeneracies occur within the pair. 

The result of this section is that the Bulk index can be expressed by means of the indices of 
the Bloch bundles of the filled pairs of bands. 

Theorem 5.7. Under the above assumption (4-3) and the gap condition ^ we have 

No/2 

l=YlliE,), (44) 



where X = T{E) is the bulk index (21) andI{Ej) is defined in Eq. (25) for the torus 



Note the bundles seen on the two sides of Eq. ( |44[ ) have different base spaces and fibers; in 
fact the dimensions of the latter are and 2, respectively. 

The proof we will give makes the simplifying assumption that the energy curves do not 
have more critical points than prescribed by the time-reversal symmetry of the Hamiltonian. 



At first, we reduce the theorem to a lemma. The base space of the bundle E in Eq. (26) can 
be extended from T to all of {{z,k) \ z e p{H{k))}. In this notation, T = T{E \ T). The torus 
¥ = 7x5"^ may then be deformed and split into A^o tori Tj, (j = 1, . . . , A'^o) each surrounding a 
pair of bands {2j — 1, 2j). Unlike T, the Tj do not need to be a Cartesian product form, since 
(43) does not imply a spectral gap uniformly in k. We then consider the bundles E^^") := E [" Tj. 



By homotopy and by the multiplicative property under splitting (Lemma 7.5), we have: 



No/2 No/2 

I{E \T)=l{[j ^(^■)) = Yl I{E^^^) , 

and Eq. (44) reduces to the following. 
Lemma 5.8. 

I{E(^^) = I{E,) . (45) 
The proof of this main lemma is deferred to Sect. [8| 



6 Quantum Hall insulators 

Much of what has been said in the previous sections has a counterpart for Hall insulators. 
In that case the results are not completely new; nevertheless they generalize the bulk-edge 
correspondence of |12j . 

The setting is the same as given by the Hamiltonians ([l| [2]) with gap condition ([s]), but 
without postulating a time-reversal symmetry. For the sake of brevity, definitions and statements 
given in the sequel rely on notations and contexts used for topological insulators. 
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Definition 6.1. [Index] Let T{ip) e GL(A^), ((/? e S^) be a continuous family. We set 



M{T) = N{Z) , 



(46) 



where T{ip) determines Z{ip) as in Eq. (16) and N{Z) is defined in Eq. (12). Alternatively, 
N{T) is the winding number of det T{ip) . 

Let E be the vector bundle with base T = x 3 {ifi,ip2) = (p and fibers of dimension 



N. On the frame bundle E{E), the right multiplication by GL{N) is defined as in (20). 



The classification of such bundles may again proceed by considering the cut torus T = 
[— vr, vr] X 5^. In fact on T there are smooth sections v : T ^ ^(^)) this is in analogy to 



Lemma |4.5[ but with simpler proof, as observed in Remark 7.1 



Definition 6.2. [Index] We set 



N{E) = N{T) 



(47) 



where T{ip2) is the transition matrix introduced in Eq. (24). See Eq. (46). The definition is 



again independent of the choice of the section v, as seen from the simplification of Lemma 4-7 
Lt should be noted that the sign of the index would flip upon interchanging + in Eq. (24). 



The index J\f{E) is just the Chern number of E, but that will not be needed. 



The definitions of bulk and edge indices parallel Defs. 5.1 and 5.2 



Definition 6.3. The bulk index is 



M = M{E) 



(48) 



where E is the bundle (26). See Eq. (4^)- The edge index Af^ is the number of signed crossings 
e of the Fermi energy /x by eigenvalues of H'^{k). They are counted positively for decreasing 
eigenvalue branches. 

Theorem 6.4. Let bulk and edge index be defined as above. Then 



Proof. The proof is contained in that of Thm. 5.4 Now its last sentence matters 



6.1 The bulk index as an index of Bloch bundles 



We consider the case of doubly periodic Hamiltonians in close analogy to Sect. 5.1 As in 



Thm. 5.7 we will express the bulk index in terms of the Bloch bundles of the filled bands. 

In the context of quantum Hall insulators it is legitimate to assume that, for fixed longitudinal 
momentum k, the bands do not overlap: 



Ai, A2, ... Ajvo > 



where Nq is again the uppermost band below the Fermi surface, cf. (44). 



Definition 6.5. The Bloch bundle Ei of the l-th band has the Brillouin zone '. 
as base and the lines [^£{k,K)] a C*^^ as fibers. 



(49) 



S^xS^ 3 {k,K) 
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Theorem 6.6. Under the above assumption (49) and the gap condition ^ we have 

N 



No 



where J\f = N'{E) is the bulk index (48) and J\f{E() is defined in (4'^) for the torus 



Proof. Consider the simphfying assumption that the energy curves do not have more critical 



points than necessary in absence of time- reversal symmetry, as explained after Eq. (50) below. 
The proof then parallels that of Thm. |5.7[ but is much simpler. In fact the complex loop is of 
the type seen in the first case of Fig. [7j but run through just once. ■ 

6.2 Bulk-edge correspondence through scattering theory 

We propose a further perspective on the bulk-edge correspondence in the doubly periodic case. It 



does not rely on decaying bulk solutions, as Thm. 6.4 did. In contrast to that result, where edge 
states are intercepted at Fermi energy, here they are right at inception, i.e. as they are born at 
band edges. That will be done by means of a result from scattering theory known as Levinson's 
theorem. In its usual form ([24j. Thm. XI.59) it computes the phase of the scattering matrix at 
thresholds. The version below computes the phase difference when a parameter is changed. 
We focus on a single band I which is separated from its neighbors. 



A^_i, > 



(50) 



and on its Bloch bundle Ei, cf. Def. 6.5, We also assume that, for fixed k e S^, the energy curve 
A = Xi{k, k) has as a function of k just two critical points, both non-degenerate, namely a max- 
imum K+{k) and a minimum K,-{k) (That assumption would not be consistent with topological 
insulators, cf. Fig. [7|) The curves K+{k) cut the Brillouin zone B into two open domains B+ 
(resp. B_) where X{k, k) is increasing (resp. decreasing) in k w.r.t. the orientation of . 



k 



t 



Figure 5: The Brillouin zone B with the domain B_ (shaded) and the curves K+{k). 



Lemma 6.7. i) There is a map r : B — > B defined by X{k, k) = X{k, r{k, n)) and the property 
that is has the two extrema k = n+{k) as its only fixed points. It interchanges and is 
real analytic in k. 
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a) There is a domain B_ 3 B_ and on there a section ^ {k,K) ^ of the Bloch bundle 
which is smooth in k and analytic in k. 

We will occasionally omit k from the notation in the sequel. In order to ensure that edge 
states are indeed "born at band edges" , we make the assumption that there are none embedded 
in the band, 

cTpp(i/«)n[A(K_),A(K+)] = 0. (51) 

The assumption is generically satisfied, but counterexamples can be constructed by taking the 
direct sum of two Hamiltonians, such that the pure point edge spectrum of one overlaps the 
band of the other. 

Lemma 6.8. At energies A = A(k), (k e B_) the edge Hamiltonian if" has a bounded eigenso- 
lution tjj'^ = {il^li{n))neN which is unique up to multiples. It satisfies 

V'i(K) =V-(K)+V+(r(K))+o(l) , (n-+c») (52) 

where the Bloch solution tp^^r^K)) 7^ is uniquely determined by '^^(/t). (We recall the relation 
(39) between ^ and Hence is a section of Ei on B_|- = r(B_). 

The bounded solution ■0" ought to be interpreted as a scattering solution for the reflection 
at the boundary n = 0. In fact, since X'{k) < on B_, represents an incoming wave, and 
■0"'" an outgoing one. 

Since B+ and B_ overlap near k+ we may introduce scattering amplitudes S+{k,K) for k 
near K.+ {k) and k > K,j^{k), resp. k < k_(A;): 

^+(r(K)) = S+[k)^- {r{K)) . (53) 

The scattering amplitudes S+{k+ + 6), which by the lemma do not vanish, play the role of 



transition matrices (24) for the line bundle E^. Therefore, by (|47j), 

M{Et)=M{S+)-M{S-) , 



where M{S+) is defined in (46). 

Definition 6.9. We say that i?" has a semi-bound state at the upper band edge A(k+) if it 
admits a bounded solution tjj'^ = {tpn)nen of H'^il)'^ = X{K+)ijjK 

Lemma 6.10. Suppose that a branch e{k) of discrete eigenvalues of H'^{k) touches the l-th band 
from above at k^, i.e. 

e{k) - X{k, K+{k)) ^ , (k^k^) . (54) 
Then H^{k:f) has a semi-bound state. 

Theorem 6.11. [Relative Levinson Theorem] Let k = ki e , {i = 1,2) not correspond to 
semi-bound states of H'^ik). Then 

Iimarg5+(A;,K+(A:) + 5) \l^= 2ttN+ , (55) 

where arg denotes a continuous argument and Nj^- is the signed number of discrete eigenvalue 
branches of H'^{k) emerging (— ) or disappearing (+) at the upper band edge, as k runs from ki 
to k2 in the orientation of . Likewise for the lower band edge and S-, except for a reversed 
count of signs in . 
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In particular the theorem may be apphed to ki = k2, i-e. to a fuh circle S^. Then it states 
J\f{S+) = N+, ((5 > 0). If ^ = A'^o is the uppermost band below the Fermi energy, then N+ =MK 
If the same assumptions hold true for all bands below it, then the bulk-edge correspondence 



is recovered in view of N_ 



N 



{i-i) 



The proofs of the results of this section are found in Sect. [91 



6.3 Comparison between two approaches 

In this section we compare the independent approaches to bulk-edge correspondence underlying 
Thms. 1131 and [6TT1 

Let H'^{k) have a semi-bound state at the isolated point A;*. Suppose that a branch e{k) of 



discrete eigenvalues disappears there into the A^o-th band from above, cf. Eq. (54) with k ] k 



whereas none emerges; see Fig. [6j By a suitable, A;-dependent energy shift we may assume that 
the upper band edge A(K+(fc)) = A-|- is constant, and pick the Fermi energy /i > A+ arbitrarily 
close to it. The branch e{k) is then decreasing and crosses the (dashed) Fermi line. For fci, /c2 
near A;* with fci < fc* < k2 we have 



Urn arg S+ {k, k+ (k) 



argdetr(A:)|^2 _ 2tt 



for ij, close enough to A+, where T is the transition matrix (36) across the Fermi line, as used in 



the proof of Thms. 5.4 and 6.4 In fact, the first equation is the specialization of Eq. (55) and 



the second one follows from the claim made after Eq. (36). In case of an emerging branch the 
r.h.s. of both equations would change sign. 





k:k 



k 



-+- 



ki k^ k2 k 



Figure 6: Left: A discrete edge eigenvalue e{k) disappearing into the uppermost band below 
the Fermi energy //. Right: Same, after energy shift. 

Both expressions on the l.h.s. arose as winding numbers. The point we wish to make here is 
that their equality can also be seen by homotopy, rather than by separate evaluation. We thus 
emphasize: 



Proposition 6.12. Under the above assumptions we have 

Iimarg5+(A:,K+(A:) + 5)|^^ = argdetr(A;)|^^ . 



(56) 
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The proof is given in Sect. [9} 



We close this section with a remark. The duahty of Sect. 6.2 reUed on the assumption that, 
for fixed k, the energy curve Xi{k,K) had a single local maximum and minimum in k. That 
naturally assigns a unique outgoing wave to a given incoming one, with transversal momenta 
n r{K). Should the assumption fail, multiple reflected waves could arise, preventing that 
essential assignment; at least without further ado, like dealing with scattering matrices instead 
of amplitudes. That failure is in fact unavoidable for time-reversal symmetric Hamiltonians (see 
Fig. [7] left). However, Prop. 



6.12 



still indicates that the index of the Bloch bundle Ej (associated 
to the pair of bands / = 2j — l,2j) can be determined by the glitches of arg S+{k, k) along the lines 
K = K+{k) of global maxima (or minima). In fact points where edge state emerge or disappear 



generically occur only on those lines, see Eq. (51 ), whence at nearby energies scattering remains 
described by a single amplitude S+. 



7 Basic bulk-edge correspondence: Proofs 



The present section contains the proofs of our basic results for topological insulators (Thm. 5.4) 



and for integer quantum Hall insulators (Thm. 6.4). We begin by supplying those relating to 
Sect. [4] (except for Lemma 4.4) and concerning the index of an abstract time-reversal invariant 
bundle. 



Proof of Lemma 4-5 The main claim (22) states that there are linearly independent vectors 



, vn{^) e such that 



1 



ii = 2,4,...,N), I.e. 



Vi{Tip) = -@Vi-l{ip) 



(57) 



We first consider either fixed point (po = ripo e T. We will prove by induction inn = 0,2,...,A^ 
that there are linearly independent vectors Vk = Vki^po) e E^^^, [k = 1, 2, . . . , n — 1, n) such 
that (57) holds true for i = 1, 2, . . . , n. Indeed, pick Vn+i such that vi, . . . , f„, are linearly 



independent and set Vn+2 



Moreover vi,. . . ,Vn+2 are still linearly independent: Xifc^i ^kVk = implies, by v 
same with A^. = X!£=i^^fc^^^! short for A = ©qA instead of A = (Ai, . . . , A„+2)- 

rank of those vectors is at least n -I- 1 we have 



Qvn+i- Then Eqs. (57) also hold for i 

im 
(Ai, 



n 



2 by 



-1. 

Qve, the 
Since the 



aA -I- oA = , 

with (a, a) e C^, (a, a) 7^ 0. Applying @q yields — aA -I- aA = and together A = A = by 
|a|2 + |a|2 ^ 0. 

The frames v{ipo) at the two fixed points may be interpolated over the line {fi = 0} x [0, vr]. 
In fact, there the bundle E is trivial (|21|. Cor. 9.5), since [0, vr] contractible to a point; and 



GL(A^) is connected. Eq. (22) is then used as a definition to extend v{<p) to a section over the 
circle {(pi = 0} x . 

The section v{(p) may then be extended further to the half torus T1/2 := [0,vr] x by 
homotopy. (We still denote by E the bundle restricted to it.) More precisely, two maps T1/2 
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Ti/2 are clearly homotopic, namely the identity map and / : {Lpi,ip2) ^ (0,(/32)- Hence F{E) 
and f*F{E) are equivalent bundles ([21], Thm. 9.4). Since the latter has fibers {f*F{E))^ = 
F{E)f(^^-j, the above section over the circle extends trivially to T1/2', by equivalence the same 
holds true for F{E). Finally, Eq. (22) extends the section v to all of T. An inspection of the 



procedure shows that the necessary smoothness of the section can be ensured. Eq. (23) then 
follows in the equivalent form eM^rip) = M{ip)e using (21). ■ 



Remark 7.1. In absence of time-reversal symmetry, see Sect.^ the existence of a global section 
on T is quite obvious. In view of the above we may just note that it exists on the circle {ipi = 
0} X S"^ = S"^ . In fact we may pick any frame t'(O) = v{2t:) and interpolate along [0, 2tt\. 



Proof of Lemma 4-6 Note that Eq. (22) states v+{^2) 
We thus have 



Qv^{—ip2)e for the boundary values. 



v+if2) = V-if2)T{ip2) = (6i;+(-(^2))er(v32) 

= {Qv-{-^2))T{-^2)eT{^2) = {e^v+{^2))eT{-ip2)eT{^2) 

Since @^ = —1 and the right action is transitive, the claim follows. 



Proof of Lemma 4.I . The two sections obey the relation v^'^^{}p) = v^^^ {}p)M{Lp) for some matrix 
M{}p) e GL(A^) continuous on T. Let M+((/92) := -/W(+7r,(^2) be its boundary values. The two 
transition matrices are then related by 

ri(992) = M_(v92)T2(v?2)M+((^2)"^ • 

More generally, we consider 

T(93i,992) := M(-V3i,992)?2(992)M((/?i,((?2)^"^ 



(58) 



as a homotopy in ^ 991 ^ vr of families in (/?2- Eq. (23) implies that, at fixed the matrices 



M(+(/9i,-) satisfy (17). Hence the family r((/?i,-) has the Kramers property. Clearly the index 



of r((/3i, •) is constant in ip\ by continuity. Note that one end of the homotopy is 

r(0,(^2) = Af(0,^2)T2((/P2)M(0,(^2)"^ , 



which by (19) has the same index as T2, while the other is (|58|). We conclude X(T2) = X(T\). 



We next come to the the proof of the main technical lemma in relation with the bulk-edge 
correspondence. 



Proof of Lemma 5.5 i). The proof of the lemma follows quite closely that of (|3], Lemma 3). 

Consider the finite difference equation (ff" — z)ip'^ = with z e C and without imposing the 
boundary condition ip^ = 0. Solutions are square-summable at n ^ +00 iff -0" = ^'"0 for 
some a e C^. Hence is eigenvalue of \I'q iff z is an eigenvalue of the operator ff", which now 
includes the boundary condition. In particular z is then real, because -ff" is self-adjoint. For 7 
as in Fig. [l]we have 7 n a{H'^) a {fi}. 

We may arrange for the absence of "flat" crossings, e(A;*) = fi, e'{k^) = 0, by adding to 
Vn{k) and Vn{k) in Eqs. (Ill N) an arbitrarily small constant. In particular, the points k^ are 
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isolated, as claimed. Moreover, they are generically simple. To show this, we perturb V, 
by tWn{k) and determine the splitting + tp, + o(i^), (t ^ 0) of a degenerate eigenvalue of 
H{kif). On general grounds Jl is found by diagonalizing W after orthogonally projecting it onto 
the unperturbed eigenspace of /i. Here, that eigenspace is the image of ker^-fi c under ^t. 
Hence the eigenvalue problem reads 

00 00 

Po[Y,i'^iTWn{K)¥^)Poa = /iPo(2(*f.)**i)Poa , (a e C^) , (59) 

n=0 n=0 



where Pq is the orthogonal projection onto ker\I/Q. By Eq. (29) the matrix in brackets on the 
r.h.s. is positive definite on C^, while that on the l.h.s. may take arbitrary Hermitian values, 
along with Wn{k^). As a result, the eigenvalues fi are generically distinct and, since e'{k^) i= 0, 
the points /c* split into non-degenerate ones. Similarly, points fc* with det ^'^ = correspond to 



/i being a Dirichlet eigenvalue for n = 1. Its perturbative splitting is determined by (59) with 



the replacement of by 1. By (29), ker^'g n ker = {0}, so that the operator on the l.h.s. can 
be chosen with independent projections under Pq and Pi. Any coincidence between eigenvalues 
of the two Dirichlet problems is thus generically lifted. 

The above argument did not pay attention to the time-reversal symmetry of the Hamiltonian 
and hence of W{}i). We do so now: A crossing 3,t — 0, vr may be displaced by perturbing 
by an arbitrarily small constant; one at /c* ^ 0,7r by letting W{k) satisfying Eq. (|4]). It still 
remains arbitrary at A;;^. ■ 

In preparation for the proof of part (ii) let us introduce the Casoratian, which is to finite 
difference equations what the Wronskian is to ordinary differential equations (both linear and 
of second order). Let H be as in Eq. (38). Given = {'ipn)nei., ^ = {^n)nez with ijjn^'fn £ 



viewed as column, resp. row vectors, let 

Cni^P,1p) = ipnA*1pn+l - <fn+lA^n ■ (60) 

Suppose that ip satisfy the Schrodinger equation in the form 

(iJ-z)V' = 0, v5(iJ-z) = 0, (61) 

where {}pll)n = V'n-i^* + V'n+i^ + '^■nVn- Then 
CI) Cn{^,ip) is independent of n e Z and denoted C{ip,ip). 



C2) Let iIj{z) solve the first Eq. (61 ). Then ip{z)* solves the second. In particular C{'tlj{z)* jtp^z)) 
is well-defined. 

C3) Items (Cl-2) apply if H,'4>,Z are replaced by ij", V^N. 
Property (CI) follows from the identity 

Cn{(p,1p) - Cn-liip,1p) = iPniAlpn-l + A*lpn+l) - {^n-lA* + ipn+lA)lpri (62) 



and Eq. (61); the others are straightforward. The Casoratian may be extended literally to ma- 



trix solutions of Eqs. (61), in which case C„,(<I>,^') is itself a matrix. Its entries are the 



24 



Casoratians of the rows and columns of resp. ^. Properties (Cl-3) hold correspondingly. 

Proof of Lemma \5.^ continued, ii) We drop k and notice that 

Co(*"(z)*,^"(z)) = *^(z)*A*^'5(z)-^'"i(z)M^f,(z) = L(z)-L(z)* . 

By (C2) the l.h.s. equals lim„^oo Cn(^''(^)*i ^"(-z)) = 0; whence the reflection property. The 
statement about the eigenvalue branch follows from part (i) and the definition of L. 

iii) We drop k = k^. Let u e be the normalized eigenvector of L{^) with eigenvalue 
= 0, whence ^'f)(//)ti = by Q. Thus, 



dz 



[u, 



dL 



/i dz 
Next we observe that 



because of (|62[ and of {H^ - z)d:,¥ = ¥. Since C„(^'''(/x)*, — > 0, (n — > (»), we 
conclude {dzl){fJ') ^ 0. Actually, equality is excluded, because Xin=i ^"('")*^"(/^) positive 
definite, as remarked earlier. 

iv) Let e{k) the eigenvalue branch crossing at A;*, whence l{e{k), k) = near fc*. The claim 
follows from 



dk 



l{e{k),k) 



5[ 
Jz 



^(k) + — 
At,fc* ok 



together with (31 ) 



This concludes the proof of Thms. 5.4 and 6.4 



7.1 Supplementary results 

This section contains a few details related to Sect. l4l some of independent interest and some 



needed in connection with Thm. 5.7, but none with the proof of the basic results, which is by 
now complete. Also addressed is the relation with other indices found in the literature. 

Polar decomposition. Given a matrix T e GL(A^), let T = PU be its (unique, left) polar 
decomposition, i.e. P = P* > 0, U*U = 1. 



Lemma 7.2. Let T e GL(A^) satisfy Eq. (15). Then so does U and the following deformations 
of T retaining that property are possible: (i) T to U , while keeping the polar part fixed; (ii) U 
to 1, while preserving unitarity; (iii) and hence T to 1. 



Proof. The right polar decomposition is T = UP with P = U ^PU. Eq. (15) states T 
@Q^T-^eo = (Gg ^C/"ieo)(G(7^P"^eo) and is hence equivalent to 



(63) 



In particular the preliminary claim holds true. 



i) Denoting by P = Ajllj, (Aj > 0) the spectral decomposition of P, the second Eq. (63) is 



equivalent to the existence of an involution i i-^ i such that = ^ and U ^IIjC/ = Gq ^IljGi 



^-1 
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By interpolating the eigenvalues but not the eigenprojections, it becomes clear that T can be 
deformed as stated. 

ii) Denoting now by C/ = Xlj ^i^i-, = 1) the spectral decomposition of U , the first Eq. (63) 
is equivalent to 11 j = 0q ^IljBo. By again interpolating eigenvalues only we arrange for U going 
to 1. ■ 



Let us also mention the following consequence: 



Remark 7.3. Let T{ip) = P{ip)U{(p) be the polar decomposition of T{ip) e GL(A^). If the 
family T has the Kramers property, then so does U. Even though the Zi{ip) are generally not the 
eigenvalues ofU{(p), it holds true that I {T) =I(U). 

Proof. By (i) of the previous lemma, T[ip) can be continuously extended to the right on an 
interval e (vr,^] in such a way that there Eq. (15) holds throughout, U{(p) = U^ir), and 
P{b) = 1. Likewise on an interval [a,0). By continuity neither index, I{U) or I{T), changes in 
the process. At this point they are manifestly equal. In fact {P | P > 0} is a cone; hence any 
continuous map P : [a,b] 3 x i-^ P{x) > with -P(a) = P{b) = 1 is homotopic to P = 1. ■ 



Proof of Lemma 4-4 ■ As mentioned, only Eq. (|18[) requires proof. Consider the families 



They enjoy the Kramers property along with Ti; they can be concatenated, since T{{tt) = T'2 (0); 
and T2 is homotopic to T2#T2 , as seen by postponing the change of M± till after that of Ti. 
Hence 

1{T2) = 1{T^#T^) = 1{T^)1{T^) , 



where I{T2) = I{Ti) by deforming M+[0) to 1 while preserving (17); and ^^{T!^) = X{M) by 
Lemma|7.2| (iii). ■ 



Alternate definition. There is an alternative way to Def. 4.8 of computing the index X{E) 
of a time-reversal invariant bundle E. Let u{(p) : T F{E) be a section satisfying (22) on the 
twice cut torus T = ([— tt, 0] u [0, tt]) x S*^, where r exchanges the left and right halves of T. The 
section gives rise to four boundary values and two transition matrices parametrized by ip2 ^ S^: 

'u(0+,(/?2) = n(0-,932)?o(v72) , n(+7r,(/72) = ^(-vr, (/?2)7V(v'2) • (64) 

Lemma 7.4. In this situation, the families Tq,Tt^ enjoy the Kramers property and 

X{E) = X{To)X{T^) . (65) 



Def. 4^ corresponds to the special case of a single cut at ipi = +7r, whence Tq = 1; a single 
cut at (^1 = could have been used there instead. 



4.6 



still applies. Let u : T — > F{E) be a 



Proof. The two families are Kramers becau se Le mma 
time-reversal invariant section as in Lemma 
v{ip) = u{ip)M{ip). It satisfies (23). Since v{ipi,p2) is continuous along ipi = 0+ and satisfies 
(|24|) along ipi = +7r we have 



and let M{ip), (ip e T) be the change of frame 



1 



M„o((^2)To(<^2)M+o(<y52) 



-1 



T{p2) = M^^{ip2)T^{^2)M+^{ip2) 



-1 



26 



where M+Q{ip2) 



M{0+,ip2), M+^(v92) = M(+7r,vP2). By Eq. ([I8j), 

1 = x(ro)x(M_oM;i) , x(r) = x(r^)x(Af_^A/-i) . 



As remarked in the proof of Lemma 4.7, the family M[—(pi, •)M((/?i, •) is Kramers for fixed 



If I e [0, vr]. It is continuous in (pi in that interval, and its index constant. Hence the claim. ■ 

Splitting the torus. We consider a time-reversal invariant bundle E on the torus T, as 
described in Sect. 4.3 We suppose moreover that (pi = +tt/2 is a distinguished pair of lines in 
the sense that 



E, 



(-f,^2) 



E, 



(66) 



(In applications, that identification of fibers may occur because the bundle E is effectively the 
pull-back of another one under a map / with /(— 7r/2, ip2) = /(vr/2, f2)-) Let us define the torus 
Ti = {{(pi,ip2) I Iv^il ^ vr/2, ip2 e S^} with identified edges pi = +vr/2, whence it still has four 
time-reversal invariant points. By ( |66[ ) Ei = E I" Ti is a well-defined time-reversal invariant 
bundle. Similarly for T2 = {{pi, P2) \ Wi — ^ ^/2, 9^2 £ 5*^} and E2 = E \ T2. 



(67) 



Lemma 7.5. [Splitting lemma] In the situation set by Eq. (gfij ) we have 

Z{E) = Z{Ei)Z{E2) . 



Proof. By Lemma 4.5 there exists a section vi{ip) of F{Ei) satisfying (22 ) on Ti, the torus Ti cut 



along ifi = +7r/2; likewise, V2{'p) on T2 with cut pi = ivr. Let Ti{{p), T2{'p) be the corresponding 
transition matrices (24). By (66) we may also introduce transition matrices M+[ip2)i {^2 £ S^) 

by 

V2{+-,IP2) = Vi{+-,(p2)M+{ip2) . 

Multiplying by G, e from the left, resp. right we obtain 

.71" X /—TT 



V2{ + -,-'P2) = Vi{ + -,-ip2)QQ M+{ip2)@0 



and thus 



eoM+(-^2) = M+{ip2)eo . 

Moreover, the l.h.s. of (68) is independent of +, whence 

Ti(v;2) = M^{ip2)M+{ip2)-^ . 



(69) 
(70) 



On T we define the section 

v{<f) 



V2{(p)M+{(p-, 



ipi\ ^ tt/2) 

ipi ^ vr) 



.1 ^ 



It is continuous by (68) and satisfies ( |22[ ): by hypothesis for \(pi\ ^ 7r/2, but also for 7r/2 ^ 
+fi ^ vr by (69). Its transition matrix is read off as 

T{ip2) = M^{ip2)T2{ip2)M+{p2y^ . 

Now Eq. ([67j) follows in the form X(r) = T{Ti)X{T2) from ([I8) [70]). ■ 
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Indices in the literature. The Z2-index has been introduced in various other forms before 



151 E] ) but to our knowledge not in the form ( 25 ) . We shah first recall a formulation [8] (there 



connected to |15j). which rests on an additional, metric structure, and later establish equivalence 
with ours when there is overlap. Like our index, it first deals with continuous families of matrices 
W{ip), (0 ^ (/3 ^ vr) which, in lieu of the Kramers property, enjoy antisymmetry = —W at 
endpoints W = W{0), W{tt). There, 

detT^ = {pfWf , 

where pf W is the Pfaffian of W; in between, consider continuous branches i'^det W{ip). One of 
them will connect pf VF(0) to I{W)piW{TT), where I{W) = +1 defines the index of the family, 
c/. ([H], Eq. (3.24)). In the special case that pf W{0) = pf W{tt), and hence det W{0) = det W{7r), 
the index reduces to 

i{w) = {-ir , (71) 

where n is the winding number of det W{ip), {0 ^ ip ^ it). 

Let us move on to time-reversal invariant bundles E. They are assumed equipped with a 
compatible hermitian metric, i.e. with an inner product (•, •) on any fiber E^p such that Q*Q = 1. 
Use is made of the fact that E is trivial if @ is disregarded. There thus is a section v{(p) of F{E) 
on the torus T 3 if = [ipi, ip2). We stress: not just on the cut torus T, but at the price of forgoing 
time- reversal symmetry. It can be taken to consist of orthonormal frames v = {vi, . . .vm)- Let 

Wij((^) = {vi{(p),Qvj{T(p)) ■ 

see dl] and ([8], Eq. (3.16) with ip, rip switched). Then W{ip)*W{ip) = 1 and W{ip)'^ = -W{Tip). 
By this last property, one may define the index as 

i{E) = i{Wo)i{w^) , 

where Wo{^p2) = W{0,ip2), W^{ip2) = W{7r,ip2), {0 ^ f2 ^ vr); see (pj, Eq. (10)). The relation 
to the indices of the present work is as follows. 

Proposition 7.6. [Relation between indices] 

i) Suppose T[(p), (0 ^ ^ vr) has the Kramers property and W{(p) := T{ip)e is antisymmetric 
at endpoints. Then 

X{T) = X{W) . (72) 
ii) For a time-reversal invariant vector bundle E as above we have 

I{E)=i{E). 

Proof. We begin with a preliminary remark. Consider the three properties of a matrix T e 
GL(iV): 

TeTe = -1 , T*T = 1 , = -eTe , 



called Kramers, unitarity and antisymmetry, the first one being indeed a restatement of Eq. ( 15 ) 
In terms olW = Te they respectively read 

WW = -1 , W*W = 1 , W^ = -W . 
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Clearly any two of them imply the third. 

i) The assumptions imply that r(0), T^n) are unitary. By Lemma 7.2 (ii) we may deform the 
family at its endpoints in such a way that T(0) = T{tt), while retaining the Kramers property 



and unitarity, and hence antisymmetry. In the process both indices (72) remain defined and 
constant. At this point, see Eqs. (13 16), I{T) = (—1)", where n is the winding number of 
detr(vp). The claim follows by ^ and detT{(p) = det W{ip). 



ii) We define a time- reversal invariant section u{(p)^ see Eq. (22), on the twice cut torus T as 



v{ip) , ("Vr ^ ^ 0) 



It has ipi = 0+ and ipi = ivr as lines of discontinuity, where Eq. (64) becomes 

by the continuity of v. Taking the inner product with the vectors of the orthonormal frames 
v{0,ip2) resp. v{iT,if2) yields Wo{ip2)£ = To{if2) and W7r(v'2)£ = 7V(v'2)- The claim follows by 
(|65l[72l). ■ 



8 The bulk index as an index of Bloch bundles: Proofs 



As a preliminary to the proof of Lemma 5.8 we consider Eq. (41) not just for |^| = 1, but for 
^ 7^ 0, together with its characteristic polynomial 

P{(, z) = det{n{0 -z) = det(^r^ +A*^ + V-z). (73) 
Its basic properties, to be proven later, are as follows. 

Lemma 8.1. P is a polynomial of degree MN in z and a Laurent polynomial in ^ of degrees N 
and —N. Moreover 

P{^ = P{t\z) (74) 
and, in the time-reversal invariant case @'H{^)@^^ = 'H{^), also 

p{tz) = Pir\z) . 

(As a function of k that case holds true for k = 0, n.) 

We then consider the Riemann surface (Bloch variety) defined by 

S = {(C,z)eC* xC|P(e,z) = 0}, (75) 

where C* = C\{0}. We make some assumptions, which are typically true: B is non-singular as 
a Riemann surface, meaning that (dP/d£^,dP/dz) 7^ at all points (^0,-^0) £ B; and, if either 
partial derivative vanishes, the corresponding second derivative does not. As a result, near a 
point [ioiZo) where dP/dz 7^ we can solve P{^,z) = as z = z(^) with z analytic; and if 
dP jd^ = 0, then 2;'(^o) = 0, but ^;"(Co) ^ 0. Similarly for interchanged roles of variables. We 
conclude that eigenvalue branches z = have only non-degenerate critical points, and only 
branch points of order 2. 
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For A; = 0, vr we will have to exceptionally allow singular points {£,o,zq) e B where both 
derivatives of P vanish. However they shall be ordinary double points, i.e. the Hessian of P is 
non-degenerate. 

Of importance is also the real Bloch variety 

Bo = {{^,z)eB\\C\ = l} . (76) 

In fact, for any z e C let m ^ be the number of with |,^| = 1 and P{(,, z) = 0, i.e. of points 
(^, z) e Bq. Then m > implies z e M, since now ^ = and P{£,, •) becomes the characteristic 
polynomial of a hermitian matrix. Moreover m is the multiplicity of z as a point of the spectrum 
(j{H). By ( |74[ ) the remaining 2N — m = 2{N — m/2) points {^^z) with |,^| ^ 1 come in pairs 
(C)-^)) {^^^^z). In particular ni is even. 

For each (^, z) e i? we consider the (geometric) eigenspace E^^^ of the eigenvalue z of the 
matrix seen in Eq. ( 73 ) . 

Lemma 8.2. 



E := {{i^z,-^) I Eq. gTJ) holds] 
is a line bundle with base B and fibers E^^z- 

In order to properly state the relation between this bundle and E we make the following 
definition. 

Definition 8.3. A generalized Bloch solution {ipn)nez of energy z and quasi-periodicity ^ 7^ 
satisfies {H — z)^ = and 

i>n+pM = i^{'4>n + pC^Tpn) , (77) 

where {iprijuez is a Bloch solution for the same ^, z. 



Eq. ( |77| ) characterizes a generalized eigenvector (of order 2) of the translation operator. In 
terms of Eq. ( 39 ) a generalized Bloch solution ■0 corresponds to ^ e C^^^ with 

-^(0$ + {A*i - AC^)C^-^ = ■ 

Lemma 8.4. Let the eigenvalue branch z = z{S^) have a critical point at ^ = and let ^{S,) be 
a local section of E at {^q^zq = z{^q)). Then ^''(^o) corresponds to a generalized Bloch solution 
for ^0,^0. 



Proof. Follows by differentiating in Eqs. (40, 41). In fact. 



(^(0-^)^' + ^* = !^ (78) 



with /(Co) = and ^dU/d^ = - AC^. 
Finally, we consider 

i?_ = {(C,z)ei?| lel <!} 

((9i?_ = Bq) equipped with 



vr : S_ ^ C, (C,z) ^ z . (79) 
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The preimage vr \z) of a point z e C consists of N points iff z e p{H). In that case we have 



where the l.h.s. is to be understood as a direct image: The fiber at z e p{H) is ©-Eg.z with 



sum over (^,2) e vr (2;); the equation itself on rests on the relation (39) between \I' e E^^^ 



and ^ e Ez- Some care applies at points (^,2;) where dP/d^ = 0, i.e. to ^ a multiple (in fact, 
double) preimage of because this seemingly results in a missing solution on the l.h.s.. It 
should however be read so as to include the there existing generalized Bloch solution. With this 
reading both sides are continuous in z e p{H). 

Likewise, the preimage 7r^^(7) of a loop 7 a p{H) is a cycle in possibly consisting of 
several loops, and we have 



■K^{E 



vr 



-1 



(7)) =E\^. 



(80) 



Proof of Lemma 5.8. We discuss a typical case only. We shall show that each side of Eq. (45) is 



associated with a cycle in or, more precisely, with a family of cycles parametrized by A; e S"^. 
Most of the discussion occurs at fixed k, which is hence once more omitted. 





(b) 






Figure 7: Left: Real energy curves Ai(e*'^), (k e S^^^l = 2j — l,2j) for various values of 
k. The first case corresponds to time-reversal invariant points k = 0,7r. Right: Loops A/(^), 
= g < 1) in the complex plane. In the first case the loop is run through twice. The second 
case is illustrated thrice (a-c) for different values of q (or k). The dashed lines represents energy 
bands of multiplicity 2, resp. 4 if thick. 
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[z,k) eJj}. For Ej 



For E^^^ = E \ Tj on the l.h.s. the association is by (80) with 7 = { 
on the r.h.s. a longer discussion is needed. The Bloch bundle Ej can be analytically continued 
from real k to k + is, as long as e is small enough. For the sake of illustration Fig. [7] shows the 
real energy curves Xi[k) = A/(e*'^), [k e S*^,/ = 2j — l,2j) for various values of k (left) next to 
the complex loops A;(^), (|^| = q) with q < I close to 1 (right). The loops are understood on 
the observation that an increasing portion of is flanked by a nearby branch running in the 
upper half-plane. As a whole, the loops may intersect the spectrum. 

In each case, the contour 7 c p{H) underlying E^^"^ may be deformed to the outermost path 
seen on the right half of Fig. [7| For concreteness, in the last case 7 is deformed into two loops, 
each surrounding one of the bands I = 2j — 1, 2j . The second loop coincides with 71 seen in 
Fig.i 



— TT 



TT 




Figure 8: Basic elements. Left: The sheet / = 2j of i3_ near Bq and the path I = /iu/2 : |^| = (7 
parametrized as ^ = e*^'^^*^), (k e S'^) for q = e^^ < 1. Right: Its image under ^ ^ ^^2^(0 is 
a cycle 71 u 72 with loops 7^ = \2j{Ii)- Further elements. Annulus A = Ai A2 = {k, \ ^ ^ 
I Im^l < e} ^ I g < 1^1 ^ 1} and its subsets Ai, {i = 1,2), where \2j is 1 to 1, resp. 2 to 1. 
The two sets are shown separated by a line (dotted). The set A2 is mapped onto the inside of 
72 ; Ai onto the complementary subset of the inside of 71 . Interpolating curve 7;^ and one of its 
preimages. Fa (both dashed); another preimage (dashed-dotted). 

We claim: 

Lemma 8.5. E \ 71 contains a subbundle E which is homotopic to the Bloch bundle Ej \ I = 
Ej. Moreover 

E = 7T^{E \T) , 

where F c 7r^^(7i) is a subcycle of the cycle 7r^-'^(7i) 1= B^. It covers 71 twice. 
As a result, F := 7r^^(7i)\F is a cycle too, and 

E := TT^{E \ f) 

provides a complementary subbundle, 

E \-ii = E®E . 
Lemma 8.6. F is contractible within B^. 
As a result E \ F, and hence E, are trivial. 
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This essentially concludes the proof of Lemma 5.8 except for those of Lemmas 8.1-8.6 In 
summary, 

^0') ^ E \ -ji^ E ^ Ej 

where = is by homotopy or by dropping trivial subbundles. One ought to worry about whether 
these equivalences persist once the dependence on fc e is included into the picture. As the 
proofs will show, the construction of L, and hence that of E, E, does not involve any choices, 
except among homotopically equivalent constructions. In particular, as k runs from to 27r, the 
cycles r match at the endpoints. ■ 

8.1 Further details 



Proof of Lemma 8.1. We expand the determinant (73) along the first row of blocks. Then, by 
the Laplace rule. 



P(^,z) = (-1 



(det A)^ 



-N 



A Vi-z A"" 




A*i 



A 



where P+ contains only powers of ^ of degree — (A^ — 1) or higher. When looking for the leading 
term in we may replace A*.^ by 0, which leaves an upper block triangular matrix. That term 
thus is (-l)*^(^+i)(det A)^ C"^- Similarly, the leading term in ^ is (-1)^(^+1) (det ^*)*^C^. 



Eq. (74) follows from det = detA^*. In the time-reversal invariant case we use detAf = 
detiQMQ-^), where 9 = diag(e, . . . , 9) in the notation of Eq. (g. We find P{^,z) = P{^, z) 
and hence the last claim. ■ 



Proof of Lemma 8.2 To be shown is that, for each {(,,z) e B, the eigenvalue z of 



cf. Eq. (41), is geometrically simple. 



If dP/dz the eigenvalue z = z{^) is even algebraically simple. In that case the eigenspace 
depends analytically on ^. If dP/dz = at (^O) -20) £ B, meaning a branch point, then we have 
by our general assumptions, 

d^P dP 

in particular z^ is an eigenvalue of 'H(Co) of algebraic multiplicity 2 and ^'(2:0) = 0, £!'[zq) ^ 0. 
It follows that nearby eigenvalues of are given as 



Z- Zo = ±s 



(81) 



, i"{zo) 

In order to show that the geometric multiplicity of zq is nevertheless 1, we denote by 11 (^) the 



projection onto the group (81) of eigenvalues z (it is analytic in ^; [16], Sect. II. 1.4), as well as 
by Ho = n(,^o) and by 

AAo = mio) - ^o)no 
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the eigenprojection and eigennilpotent of zq. The claim amounts to A/q 0, and we prove it by 
contradiction: A/q = imphes that 

? - ?0 



has a removable singularity at ^ = Co- Hence z — zq = 0{(^ — ^q), in violation of (81 ) 



The proofs of Lemmas 8.5 and |8.6| will mostly deal with a single at a time. We will give 



them first for the last case of Fig. [7] and adapt them later to the others. 



Proof of Lemma 8.5, first part. We consider a homotopy of curves 7a, (2 ^ A ^ 1) which lie 
between 72 and 71. In particular, we keep the endpoints of 7a fixed, as indicated by the dashed 
line in Fig. [s] (right). Among the many curves in t^~^{'^\) let us track the curve, F^ c 
arising from /2 at A = 2 by continuity, i.e. F2 = /2- The curve Fa has to stay outside of A. 
Indeed, it can cross neither I2, as T^iJ'x) would end up inside 72, nor /i, as ■7r(FA) would approach 
71 from the outside. As Fa does not penetrate into A, it does not approach the real k axis {i.e. 
1^1 = 1) on the branch of A2j. Moreover it does not either on any other branch of Bq, since those 
map to other, disjoint bands. That ensures that Ii u Fa is a cycle in B- and that the bundle 
E \ (/iuFa) is defined and continuous in A. Finally, the homotopy of bundles -K^i^E \ (/i uFa)), 
(2 ^ A ^ 1) does the job: For A = 2 we have 

7r,(^ \{hvjh))=Ej \I 

by construction; for A = 1, 7r(/i u Fi) = 71 u 71 = 71 and 

7r,(^ \ [h u Fi)) d n,{E \ 7r-\-fi)) = E \ 



by (80). Thus Lemma 8.5 holds with F = /i u Fi. 



Let us recall the count of preimages 7r^^(2;) done in relation with Eq. (79). We find it conve- 
nient to pretend that all points z e C, i.e. including z e (t{H), have precisely preimages under 
TT. This is enforced by including the m preimages on Bq, but by counting them with weight 1/2, 



as explained in the sequel of Eq. ( 76 ) . The preimages of a loop intersecting the spectrum then 
is a pseudo-cycle, where it is tolerated that curves break up at pairs of points in Bq with same 
Xi{k) (see the joined dashed and dashed-dotted curves in Fig. Isl left). 



Proof of Lemma \8.6[ We first observe from the previous proof that any point z inside 71 has 
N — 2 preimages z) outside of and away from Ii u Fi. We then extend 7a from 1 ^ A ^ 2 
to 2 ^ A ^ 3 so that it contracts 72 to a point 73 = {zq} inside of it. We then consider the 
deformation Fa, (1 ^ A ^ 3) of pseudo-cycles covering 7a twice and starting with Fi = /i u Fi. 
By continuity in A, Fa moves inside of Ii u Fi, while N — 2 preimages of 7a remain outside. As 
a result Fa := vr^-'^(7A)\FA stays away from Bq and is a cycle in ■ 



Proof of Lemma 8.5. secondpart. We extend the proof to general values of /c. They are illustrated 
by the three cases seen in Fig. [T] with the middle one having three variants (a-c). The argument 
will rely on two sheets {I = 2j — 1, 2j) of as opposed to just one before (see Fig. [8] left). The 
sheets share two branch points: In the first case (A; = or tt) they rather are, to be precise, two 
double points at k = and vr. In the variant (a) each of them gives way to two branch points. 
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out of which one is in and actually in the annulus A; in (b) the branch points have reached 
the path / bounding A (this might occur at different k for the two points), and in (c) they have 
left A. This is seen in Fig. [7] (right) as a rearrangement of loops through an intersection point 
at (b). 

In the first case, where the curves do not intersect the spectrum, the conclusion is reached as 
before, but without need for deformations. The last case has been dealt with before. We thus 
turn to the three variants (a-c) of the middle case. 

(a) We deform the inner parts 72 of the curve till they reach the outermost parts 71, while 
keeping the intersection points fixed. By the same arguments as before the preimage selected 
by continuity remains outside the annuli A of both sheets. 

(b) The cycles in (a) and (c) are related by the insertion of a loop (in z) surrounding the 
intersection point. Its preimage on the two sheets (in ^ or «;) surrounds the branch point twice. 



see Eq. (81 ). After one turn the vectors ^2j-i(^i 1^) and "^2j{k, k) get exchanged along with the 



eigenvalues z. Each of them defines in two turns what their span [^2j~iik, k), ^2jik, k)] defines 



in just one. That 2-dimensional fiber has a limit at the branch point, as remarked after Eq. (81), 
and unlike the two vectors individually. As a result the insertion is trivial. 

(c) The two sheets of the pair give rise to two curves, which are mutually and self-intersecting. 
First, let us pick one curve and expand its inner loop till it reaches the outer loop or the other 
curve. Till then the preimage of that loop avoids the set A of its own sheet, but will reach from 
the outside that of the other sheet. (This is seen from the previous item.) At this stage the 
parts crossing the spectrum are run through twice but oppositely in B, and may be removed. 
We proceed similarly with the other curve. The resulting cycle in B no longer projects on the 
spectrum under vr. 

Finally we observe that the deformations in (a) and (c) are smoothly connected to one an- 
other by the move explained in (b), but also with the first and last case respectively. The 
deformations are thus continuous function of A;. ■ 



Likewise the proof of Lemma 8.6 extends to the other cases, too 



9 Bulk-edge correspondence through scattering theory: Proofs 



Proof of Lemma \ 6.T\ i) The map r(«;) = r{k, k) is well-defined by the required properties, 
and it remains to show that it is real analytic in k. Away from critical points, k k+, let 
F{K,r) = — A(r), so that F{K,,r{K)) = 0. Since dF/dr = — A'(r) 7^ for r = ^(k), the 
claimed analyticity follows by the implicit function theorem. Near a critical point, say k+, the 
argument must be modified. Let there 



A'(ai), (k = r). 



Note that F(/t,r) = still has the solution r = r(«;), but no longer r = k, except for k = k+. 
Since 

dF _ X{k) - A(r) - A'(r)(K - r) 
dr (k — r)-^ 

equals A"(k_|_)/2 7^ at (k, r) = (K+,r(K+)) the solution r = r(«;) is analytic also near k = k+. 
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ii) We recall the Def. 6.5 of the Bloch bundle E^. As shown in Fig. [sj ]B_ does not contain any 
non-contractible loop winding around k e S^. Thus there is no obstruction for a smooth section 
^{k, k) ^ of E£ on IB-, or even on small complex neighborhood in k thereof. To be shown is 
that there is one, ^^{k, k), which at fixed k is analytic in k for (/c, k) in that neighborhood. As 
a preliminary, let P{k,K) be the Riesz projection onto the fiber {E()k^^ e C^*^, 



1 



where = 'H{e'''^) is defined in Eq. (41), the variable k is suppressed, and the integration 

contour surrounds once the eigenvalue Xi{k) of T-1{k), and no further ones. We observe that P 
is analytic in k. 

We construct first by setting ^^{k, k.) = ^{k, k) for, say, k = fi;-|-(/c); then by extending 
it in K (at fixed k) through parallel transport: 



Pdn^- = , 
This is feasible. In fact, the ansatz 

^-(k) = ^(A;)A(k) 

reduces both equations to 



(A(M e 



(82) 



p{d^^) + {p^)d,\ogx = o . 

Here = /(P^*) for some function /(k) e C, since the fibers are lines. The equation 

thus reads d^logX = — /, which can be integrated starting from A(k+) = 1. It remains to 
verify the Cauchy-Riemann condition d^^^ = 0. By the Eqs. (82) and d^P = we indeed have 

e;^- = p(d;^-) = o. ■ 



In preparation for the proof of Lemma 



6.8, or actually of an extension thereof, let us list two 



further properties of the Casoratian, Eq. (60). This time, the energy z in Eq. (61) is real. 



C4) Let he a, Bloch solution of quasi-periodicity ^ = e*'', (k e S^) and energy z = A(k). Then 

M-l 

C(V'*,V') = -iA'(K) 2 CV'n. (83) 

n=0 

C5) Let ip = iP{k) be a Bloch solution as in (C4) for k near k+. Then ip'[K,+ ) is a generalized 



Bloch solution by Lemma 8.4 and 



M-l 



(84) 



n=0 



The derivation is as follows. We recall that ^ in Eq. (39) satisfies Eq. (41) with 'H(^) = 'H(^)'' 
Eq. (78) thus yields the Feynman-Hellmann equation 
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by ^{d/d£) = —id/dK, the inner product being the standard one on . The l.h.s. equals 



as claimed in (C4). Differentiation of ([83j) yields (|84j) because of C{4)*,(p*) = -C{(p,4))* and 
X'{k+) = 0. 

Let us consider a complex neighborhood of the band A under consideration: 

N{6o) = {X{K)\\lmK\<6o} , {60 > 0) 

The following lemma manifestly implies Lemma |6.8[ which therefore will not require a separate 
proof. 

Lemma 9.1. [Solutions of at most small exponential growth] For small enough Eq > there is 
5o > such that we have: For z = A(/t) e N{do) there is (up to multiples) a unique solution 

V'M 0/ 

{H^ - z)i^^ = , i>l = o 
with ipt = 0{e^°"'), (n +00). It is of the form 

4 = V'i^') + 0{e-'n , (n - +«) (85) 



where 
and 



i^l^'^ = i^^^ + ^1^) ^ (86) 



i) for z ^ X[k+): ipn \ tpn^ ire Bloch solutions for H with = e*'^, resp. ^ = e*''^"-'; 

a) for z = A(k+); Vn^^ V'n^^ '^re Bloch, resp. generalized Bloch solutions with ^ = e^'^ = e^^^'^\ 

Moreover, if z = X{k), {k e S^), then ipn = 0(1) C'^d tpn\ipn^ 9^ in case (i); and ipn = 0{n) 
in case (ii). 

Proof. At first, let z belong to the band under consideration, i.e. z = X{k), (k e 5"*^). As 
explained in connection with the Bloch variety Eq. (75), we have P{S,,z) = for m = 2 values 
of with |,^| = 1 (counting multiplicities) and for A^ — 1 values with |^| ^ e^^" and small enough 
eo > 0. By possibly making it smaller the same remains true for z e N{6o) and small 5o > 0, 
provided |^| = 1 is replaced by e^^° < |^| < e^°. Any eigensolution oi H, which remains 
0(e^°"), (n — > +00) is thus of the form 

N+l 

V^„ = Vl')+Vi')+ 2 (rieZ) (87) 

i=3 

where (ipn^) is a (generalized) Bloch solution corresponding to the above ^'s, in the stated order. 
Disregarding the boundary condition, the same applies to eigensolutions "0" of -ff" for n > uq, 
by means of the map ( 28 ) . 

We claim that up to multiples there is precisely one such solution ■0", once the boundary 
condition -^q = is imposed. Without any conditions the space of solutions V'" has dimension 
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2N . The subspaces determined by 'i/'o = 0' resp. 0(6^""), have dimensions A^, + 1. Their 
intersection V thus has dimension at least 1. We make two claims: (a) ynj^/'" | if)^^'^'^ = 0} = {0}, 



cf. (85), which implies Eq. (86) and dimT/ ^ 2; and the stronger: (6) diml/ ^ 1, and hence 
dimy = 1, which is the statement of uniqueness of V'"- Properties (a, 6) are stable, because the 
nullity dim(Vi n V2) of a pair of subspaces is an upper semi-continuous function of them (|16j. 
Thm. IV. 4. 24); hence it suffices to prove them for z = A(k) with k real. 

We begin with (a): The opposite would amount to an embedded eigenvalue, which is ruled 



out by (51 ). To prove (6), together with the remaining claims, we consider first the case k ^ k+, 
where ^ = e*'^,e*^('^) are different. Suppose, indirectly, dimy = 2. Then, by taking a suitable 
linear combination of solutions we could arrange for V^^^ ^ 0, V'^^) = (or viceversa) . This 
would imply the contradiction 



= Cn 



lim C„('0«*,'0") = C(V'(i)*,V'^^)) ^ 



by V'o ^ ^^"^ preliminary remark (C4). Finally we are left with k = (or This time 
dimV = 2 would imply that there exist two solutions V'a, '0^ ^ o'^^ with ?/;a = "0) V'a = 



in Eq. (87) and the other with ■0^ 



(1) 



0, 



. The resulting contradiction is 



o = Co(V'!,*,V'l) = C(V''*,V^)^o 



by (C5) and A"(k+) 7^ 0. Moreover, the bounds 0(1) and 0{n) follow from Eqs. (40, 77) with 
lei = \e'^\ = 1. ■ 



Proof of Lemma 6.10. Let again k be fixed till further notice. We consider z in a complex 



neighborhood of the upper band edge A(k+). There the expanded part of the solution (85) may 
alternatively be written as a linear combination of 



. 5k ' 



(k 7^ K+J 



where iP~{k) is the section of Bloch solutions of Lemma 6.7 (ii), rather than as linear combination 
of Bloch solutions The advantage of the basis {^^'^\k),iIj^'^\k)] is that it does 

not degenerate as k ^ k+. The solution (85) is unique up to a multiple, which we shall fix by 



means of a prescription independent of k. For instance, since 7^ 0, there is a linear functional 
£ on such that ^ for K near k^. We impose the normalization i{ijj\{K)) = 1. We 

then have 

V'I(k) = a(K)Vi') (/^) + (k) + 0{e-'n , (n ^ +(») (88) 

with a, /3 analytic near k^. Moreover, a,/3 are even under r since are. A semi- 

bound state is tantamount to /3(k_|-) = 0. For k 7^ we may also write 



where 



4(k) = /(aj)^, 

/(k) = a(K) -h 



K) + f{r{K))^l:-{T{'^)) + 0{e-'n 
/3(k) 



r(K) 



(89) 
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is analytic in a punctured neighborhood of k+. Comparison with Eqs. (52, 53) yields for later 
use 

firiK)) 



S+(k) 



[K> K: 



(90) 



- Im K 



For luiK < the Bloch solution ipni'^) exponentially diverging for n — > oo due to |^| = e" 
Based on z = X{k), edge state energies e > A(k+) close to A(k+) occur iff /(k) = for some k 
close to with ImK < 0. For that to happen /3[k) has to be correspondingly small: 



/3(k) = a(K)(r(K) — k) 



(91) 



Let us now reintroduce the parameter k of the lemma and of its assumption Eq. (54). Then 
/3(A;*, K+(fc=f)) = follows by continuity. This concludes the proof. ■ 



Proof of Theorem 6.11. Let us first deal with the simple case that the interval [ki , /C2] does not 

and 



contain any /c^ where H'^{k<^) has a semi-bound state. Then A^_|- vanishes by Lemma 
so does the l.h.s. of Eq. ([55|). In fact, by 



6.10 



90) we have 



liin 5(A;, K,+ {k) + 5) 



-1 



uniformly in A; e [A;i,/i:2], since l3[k, K+^k)) ^ 0. We consider next the case where such points 
/c* are present. Though they are generically isolated, we will treat the general case, where they 
form intervals / c [^1,^2]. Those are closed, countably many and possibly consisting of single 
points. By the first case and by compactness it suffices to consider such an interval / and ki, k2 
sufficiently close to its endpoints. 

Let first ke I. Then a{k,K+{k)) 7^ 0, (/c e I) since (3{k,K+{k)) = 0; see Eq. (88). We claim 



there is a punctured disk of fixed radius centered at k_|-(A;), which remains free of zeros k of 



f{k,K). Indeed, if that radius is small, Eq. (91) can not hold true there, because its r.h.s. is 
linearly large in n — K-^.{k) 0, while the l.h.s. is quadratically small, as /3 is even. The claim 
extends by continuity to A: e [ki , k2] , provided the punctured disk is replaced by an annulus 
A{k) of fixed radii. 



We recall that by (51) real k 7^ n+ik) do not occur as zeros either. Let N(k) be the number 
of zeros inside the annulus and having Iuik < 0. Thus 

-~N+ = N{k2) - N{ki) . 



Eq. (55) then follows from the claim that 



YiuisxgS+{ki,K+{k,) + 5) = -2^{N{ki) + -) , (i = l,2) 



(92) 



In proving it we drop fcj from the notation, e.g. A = A{ki). We note that by assumption 
/3(k+) 7^ 0, whence there is a disk D centered at k+ which is free of zeros of /. In particular, 
K-f- + 6 e D for small 6 > 0, and k+ + A e ^4 for suitable A > likewise for the images under 
r, see Fig. [9j 

We consider the contour A u F u r(A) u 7, where A = + 5, k+ + A], F 1^ ^ n {Im k < 0} 
joins K+ + A to r{K-\- + A), and j ^ D n {Iuik < 0} joins r{K+ + 5) to k+ + 6. Let us denote 
by C that contour with opposite (positive) orientation. By the argument principle we have 



fi^r'^/d^ 

OK 
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Figure 9: The disk D, the annulus A and the contour C in the K-plane. 



The contour may also be spht into two parts, and their contributions computed otherwise. First, 



by (90) 



aig S{k+ + 6) = arg /(r(K+ + 6)) — arg + (^) = Im 

J/ 



OK 



Aurur(A) 

(Note that the first equahty is consistent with arg being a continuous argument in /c e [ki , ^2] ; 
so is the second, because the path remains free of zeros as k changes, unhke 7). Second, by 



Eq. (89) 



hm {(k) 



-1^ 



dn = in . 



Together, this proves (92) and hence the theorem for the upper band edge. The case of the 



lower band edge is similar, except that incoming states are found at k < K^{k). This explains 
the reversed count of signs in . ■ 



Proof of Proposition 6.12. We write k only when necessary. Let us recall that z = A(k) 



bijectively maps a neighborhood of /t+ in the half-plane Im k < to one of the band edge A-|- in 
the (complex) complement of the band. It will be understood that z and k are so related. We 
first consider the case of a disappearing branch and reformulate the statement using the notation 



from the proofs of Thms. 6.11 and 5.4 Since no eigenvalue branch is present at ^2 > A;*, the 



two sides of Eq. (56) are, after dividing by 27r, those of 



r/(A;i,K)-ld/ = -^ f l{z,k)-^dl, 
2m Jc 2m Jg^ 



(93) 



where 



C is the contour in the K-plane described in Fig. |9j 

dD 1= T is the contour encircling the crossing point with the Fermi line and described in 
Fig. [4} 



l{z,k) is the eigenvalue of L{z,k) described in Lemma 5.5 (ii). (We recall that the use 



made of L{z,k) rested on Eq. (30)) 



L{z) has an eigenvalue l{z) = iff /(k) = (with both zeros being of first order), since both 



conditions are equivalent to 2; = e. Hence Eq. (93) holds provided the contours are homotopic 
under z = X{k). 
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To show this, let us visuahze dD cz T in the setting of Fig. [T] That contour can be rotated 
without intersecting the discrete eigenvalue branch, and so as to lie in a 2;-plane at fixed k = 
ki. The contour remains positively oriented w.r.t. the orientation of that plane and, after 
substitution z = X{k), it is homotopic to C at /ci < A;,^. 

In the case of an emerging eigenvalue branch, the l.h.s. of (93) is evaluated at k2 > k^ 



and acquires a minus sign, cf. (56). However, when dD is rotated as just prescribed it ends up 
negatively oriented in the z-plane at k = k2- Hence the modified Eq. (93) still holds. ■ 
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